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ABSTRACT 

Gas and star velocity dispersions have been derived for eight circumnuclear star- 
forming regions (CNSFRs) and the nucleus of the spiral galaxy NGC 3310 using high 
resolution spectroscopy in the blue and far red. Stellar velocity dispersions have been 
obtained from the Call triplet in the near-IR, using cross-correlation techniques, while 
gas velocity dispersions have been measured by Gaussian fits to the H/3A4861 A and 
[Om] A 5007 A emission lines. 

The CNSFRs stellar velocity dispersions range from 31 to 73kms _1 . These values, 
together with the sizes measured on archival HST images, yield upper limits to the 
dynamical masses for the individual star clusters between 1.8 and 7.1 x 1O 6 M , for 
the whole CNSFR between 2 x 10 7 and 1.4 x 10 8 M , and 5.3 x 10 7 M for the nucleus 
inside the inner 14.2 pc. The masses of the ionizing stellar population responsible for 
the Hll region gaseous emission have been derived from their published Ha luminosities 
and are found to be between 8.7 x 10 5 and 2.1 x 10 6 M for the star-forming regions, 
and 2.1 x 10 5 M for the galaxy nucleus; they therefore constitute between 1 and 7 
per cent of the total dynamical mass. 

The ionized gas kinematics is complex; two different kinematical components seem 
to be present as evidenced by different line widths and Doppler shifts. 

Key words: HII regions - galaxies: individual: NGC 3310 - galaxies: kinematics and 
dynamics - galaxies: starburst - galaxies: star clusters. 



1 INTRODUCTION 

The gas flows in disc of spiral galaxies can be strongly per- 
turbed by the presence of bars, although the total disc star 
formation rates (SFR) do not appear to be significantly af- 
fected by them (Kennicutt 1998a). These perturbations of 
the gas flow trigger nuclear star formation in the bulges 
of some barred spiral galaxies. External environmental in- 
fluences can have strong effects on the SFR, among them, 
the most important by far, is tidal interactions. The young 
extragalactic star clusters belonging to these systems have 
been the aim of different studies during the last decades 
(e.g. Diaz et al. 1991; Whitmore et al. 1993; Mengel et al. 
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2005; Bastian et al. 2005, 2006; Mengel et al. 2008). The en- 
hancement of the SFR is highly variable depending on the 
star formation conditions, the degree of enhancement rang- 
ing from zero in gas-poor galaxies to around 10-100 times 
in extreme cases (Kennicutt 1998a). Much larger enhance- 
ments are often seen in the circumnuclear regions of strongly 
interacting and merging systems (Kennicutt 1998a, b). 

Yet, the bulges of some nearby, non-interacting, spi- 
ral galaxies show intense star-forming regions located in a 
roughly annular pattern around their nuclei. In the middle 
of last century, Morgan (1958) classified a sample of galax- 
ies using as the principal classification criterion the degree 
of central concentration of light of each galaxy. An appar- 
ent fairly common phenomenon in some types of galaxies 
was pointed out by Morgan: their nuclear regions can con- 
sist of an extremely bright, small nucleus superposed on a 
considerably fainter background, or it may be made up of 
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multiple "hot-spots". Almost a decade later, Sersic & Pas- 
toriza (1965) suggested a relationship between the existence 
of a bar and the presence of abnormal features in their nu- 
clei for a survey of bright southern galaxies. These authors 
extended the survey to the whole sky (Sersic & Pastoriza 
1967), and found that ~14% of these galaxies presented 
peculiar nuclei. The distinctive nature (with respect to the 
more extended star formation in discs) of the luminous nu- 
clear star-forming regions was fully revealed with the open- 
ing of the mid- and far-infrared (IR) spectral ranges (see e.g. 
Rieke & Low 1972; Harper & Low 1973; Rieke & Lebofsky 
1978; Telesco & Harper 1980). 

In general, CNSFRs and giant Hn regions in the discs of 
galaxies are very much alike, although the former look more 
compact and show higher peak surface brightness (Kenni- 
cutt et al. 1989) than the latter. CNSFRs, with sizes going 
from a few tens to a few hundreds of parsecs (e.g. Diaz et al. 
2000) seem to be made of several Hll regions ionized by lu- 
minous compact stellar clusters whose sizes, as measured 
from high spatial resolution Hubble Space Telescope (HST) 
images, are seen to be of only a few parsecs. Their large 
Ha luminosities, typically higher than 10 39 ergs _1 , point to 
relatively massive star clusters as their ionization source. 
Although these Hll regions are very luminous (M v between 
-12 and -17) not much is known about their kinematics or 
dynamics for both the ionized gas and the stars. It could 
be said that the worst known parameter of these ionizing 
clusters is their mass. As derived with the use of population 
synthesis models their masses suggest that these clusters 
are gravitationally bound and that they might evolve into 
globular cluster configurations (Maoz et al. 1996). Further, 
deeper analysis as to whether or not such cluster 
would survive the hostile environment of the cir- 
cumnuclear regions is extremely interesting but lies 
outside the scope of the present work. Classically it is 
assumed that the system is virialized hence the total mass 
inside a radius can be determined by applying the virial the- 
orem to the observed velocity dispersion of the stars (a*). 
As pointed out by several authors (e.g. Ho & Filippenko 
1996a), at near-IR wavelengths (8500 A) the contamination 
due to nebular lines is much smaller and since red supergiant 
stars, if present, dominate the light where the Call AA8498, 
8542, 8662 A triplet (CaT) lines are found, these should be 
easily observable allowing the determination of a* (Terlevich 
et al. 1990; Prada et al. 1994). 

The equivalent width of the emission lines are lower 
than those shown by the disc Hll regions (see e.g. Kennicutt 
et al. 1989; Bresolin & Kennicutt 1997; Bresolin et al. 1999). 
Combining GEMINI data and a grid of photo-ionization 
models Dors et al. (2008) conclude that the contamination 
of the continua of CNSFRs by underlying contributions from 
both old bulge stars and stars formed in the ring in previ- 
ous episodes of star formation (10-20 Myr) yield the observed 
low equivalent widths. 

This is the third paper of a series to study the pecu- 
liar conditions of star formation in circumnuclear regions 
of early type spiral galaxies, in particular the kinematics 
of the connected stars and gas. In this paper we present 
high-resolution far-red spectra and stellar velocity dispersion 
measurements (a*) along the line of sight for eight CNSFRs 
and the nucleus of the spiral galaxy NGC3310. 

NGC3310 (UGC5786, Arp217) is a starburst galaxy 



classified as an SAB (r) be by de Vaucouleurs et al. 
(1991), with an inclination of the galactic disc of 
about i^40 (Sanchez-Portal et al. 2000). Its coordinates 
are a 2 ooo = 10 h 38 m 45!9, £ 20 oo=+53° 30' 12"(de Vaucouleurs 
et al. 1991). These authors derived a distance to the galaxy 
equal to 15Mpc, giving a linear scale of ~ 73 pc arcsec -1 . 
This galaxy is a good example of an overall low metal- 
licity galaxy, with a high rate of star formation and very 
blue colours. This galaxy has a ring of star forming regions 
whose diameter ranges from 8 "to 12 "and shows two tightly 
wound spiral arms (Elmegreen et al. 2002; van der Kruit 
& de Bruyn 1976) filled with giant Hll regions. These cir- 
cumnuclear regions present low metal abundance (0.2-0.4 Zq 
Pastoriza et al. 1993), in contrast to what is generally found 
in this type of objects, which show high metallicities (Diaz 
et al. 2007). In fact, in most cases, the [Om] A 5007 A for- 
bidden line can barely be seen (see e.g. Hagele et al. 2007, 
2009, hereafter Paper I and Paper II, respectively). 

The ages indicated by the colours and magnitudes of 
the star formation regions are lower than 10 Myr (Elmegreen 
et al. 2002). From near-IR J and K photometry these au- 
thors derived an average age of ~10 7 yr for the large scale 
"hot-spots" (star forming complexes). From the observed 
CaT line in the Jumbo Hll region Terlevich et al. (1990) 
derived an age around 5 to 6 Myr. Elmegreen et al. (2002), 
comparing their data with Starburst99 models (Leitherer 
et al. 1999), estimated masses of the large "hot-spots" rang- 
ing from 10 4 to several times 1O 5 M0. They found 17 can- 
didate super star clusters (SSCs) with absolute magnitudes 
between Mb =-11 and -15 mag, and with colours similar to 
those measured for SSCs in other galaxies (see for example 
Barth et al. 1995; Larsen et al. 2001). 

We have measured the ionized gas velocity dispersions 
(a g ) from high-resolution blue spectra using Balmer H/3 and 
[Om] emission lines. The comparison between cr* and a g on 
an ample sample of objects might throw some light on 
the yet unsolved issue about the validity of the gravitational 
hypothesis for the origin of the supersonic motions observed 
in the ionized gas in Giant Hll regions (Melnick, Tenorio- 
Tagle & Terlevich 1999). In Section 2 we describe the obser- 
vations and data reduction. We present the results in Section 
3, the dynamical mass derivation in Section 4 and the ion- 
izing star cluster properties in Section 5. We discuss all our 
results in Section 6. Finally, the summary and conclusions 
are given in Section 7. 



2 OBSERVATIONS AND DATA REDUCTION 

The data were acquired in February 2000 using the two 
arms of the Intermediate dispersion Spectrograph and Imag- 
ing System (ISIS) attached to the 4.2-m William Herschel 
Telescope (WHT) of the Isaac Newton Group (ING) at the 
Roque de los Muchachos Observatory on the Spanish island 
of La Palma. The CCD detectors EEV12 and TEK4 were 
used for the blue and red arms with a factor of 2 binning in 
both the "x" and "y" directions with resultant spatial reso- 
lutions of 0.38 and 0.36 arcsec pixel -1 for the blue and red 
configurations respectively. The H2400B and R1200R grat- 
ings were used to cover the wavelength ranges from 4779 to 
5199 A (A c = 4989 A) in the blue and from 8363 to 8763 A 
(A c = 8563 A) in the red with resultant spectral dispersions 
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Table 1. Journal of Observations 
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of 0.21 and 0.39 A per pixel respectively, providing a com- 
parable velocity resolution of about 13kms -1 . A slit width 
of 1 arcsec was used which, combined with the spectral dis- 
persions, yielded spectral resolutions of about 0.4 and 0.7 A 
FWHM in the blue and the red, respectively, measured on 
the sky lines. Table 1 summarizes the instrumental configu- 
ration and observation details. 

The data were processed and analyzed using iraf 1 rou- 
tines in the usual manner. Further details concerning each 
step can be found in Paper I. With the purpose of measur- 
ing radial velocities and velocity dispersions, spectra of 11 
template velocity stars were acquired to provide good stellar 
reference frames in the same system as the galaxy spectra 
for the kinematic analysis in the far-red. They correspond to 
late-type giant and supergiant stars which have strong CaT 
features (see Diaz et al. 1989). The spectral types, luminos- 
ity classes and dates of observation of the stellar reference 
frames used as templates are listed in Table 2 of Paper I. 



3 RESULTS 

Two different slit positions (SI and S2) were chosen in order 
to observe 8 CNSFRs and the nucleus of the galaxy. One of 
them the conspicuous Jumbo region, labelled J, is the same 
region labelled R19 by Diaz et al. (2000) and region A of 
Pastoriza et al. (1993). Balick & Heckman (1981) dubbed 
it Jumbo, given its extreme luminosity (100 times more lu- 
minous in IR than 30 Dor in the Large Magellanic Cloud, 
Telesco & Gatley 1984). 

Fig. 1 shows the selected slits, superimposed on photo- 
metrically calibrated optical and IR images of the circumnu- 
clear region of this galaxy acquired with the Wide Field and 
Planetary Camera 2 (WFPC2; PCI) and the Near-Infrared 
Camera and Multi-Object Spectrometer (NICMOS) Camera 
2 (NIC2) on board the HST. These images have been down- 
loaded from the Multimission Archive at STScI (MAST) 2 . 
The optical image was obtained through the F606W (wide 
V) filter, and the near-IR one, through the F160W (H). The 
IR image does not cover the whole circumnuclear ring, and 
it does not include the regions Rl, R2, R7 and X. The CNS- 
FRs have been labelled following the same nomenclature as 
in Diaz et al. (2000), with the nomenclature given by Pas- 
toriza et al. (1993) for the regions in common [E for R4 
and A for R19 (the main knot of the Jumbo region)] within 

1 IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Sci- 
ence Foundation (NSF). 

2 http:/ /archive. stsci.edu/hst/wfpc2 
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Figure 2. F658N (narrow-band [Nil] filter, at the redshift of 
NGC 3310, z = 0.0033 (Haynes et al. 1998), equivalent to the Ha 
narrow band filter) image centred on the galaxy obtained with 
the ACS camera of the HST. The orientation is north up, east to 
the left. The location and P. A. of the WHT-ISIS slit positions, to- 
gether with identifications of the CNSFRs extracted, are marked. 
The radii of the circles, centred at the position of the nucleus, are 
8 and 18 ". 

parentheses. In both studies the regions observed are identi- 
fied on the Ha maps. We have also downloaded the F658N 
narrow band image (equivalent to Ha filter at the redshift 
of NGC 3310) taken with the HST Advanced Camera for 
Surveys (ACS), shown in Fig. 2. The plus symbols in both 
panels of Fig. 1 and in Fig. 2 represent the position of the 
nucleus as given by Falco et al. (1999). 

Fig. 3 shows the spatial profiles in the H/3 and 
[Oiii] 5007 A emission lines (upper and middle panels) and 
the far-red continuum (lower panel) along each slit position. 
Due to the presence of intense regions (J in the blue and the 
galaxy nucleus in the red range in the case of SI, and R5+R4 
in the blue range for S2) the profile details are very difficult 
to appreciate, therefore we show some enlargements of these 
profiles in Fig. 4. In all cases, the emission line profiles have 
been generated by collapsing 1 1 pixels of the spectra in the 
direction of the resolution at the central position of the lines 
in the rest frame, A 4861 and A5007A respectively, and are 
plotted as dashed lines. Continuum profiles were generated 
by collapsing 11 resolution pixels centred at 11 A to the blue 
of each emission line and are plotted as dash-dotted lines. 
The difference between the two, shown by a solid line, cor- 
responds to the pure emission. The far-red continuum has 
been generated by collapsing 11 pixels centered at A 8620 A. 
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Figure 1. Left: F606W (wide V) image centred on NGC 3310 obtained with the WFPC2 camera (PCI) of the HST. Right: HST-NICMOS 
(NIC2) image obtained through the F160W filter. For both images the orientation is north up, east to the left. The location and P.A. of 
the WHT-ISIS slit positions, together with identifications of the CNSFRs extracted, are marked. 



The regions of the frames to be extracted into one- 
dimensional spectra corresponding to each of the identified 
CNSFRs, were selected on the continuum emission profiles 
both in the blue and in the red. These regions are marked 
by horizontal lines and labelled in the corresponding Fig- 
ures. In the H/3 profiles we find two almost pure emission 
knots, one for each slit position, labelled Y and X by us (see 
Fig. 2), respectively. The former of these regions seems to 
be located at the tip of one vertical arm formed by pure 
emission regions, since it can be easily appreciable in the 
Ha image from the ACS (Fig. 2) but it is almost invisible 
in the WFPC2-PC1 V-band image (see Fig. 1). 

The spectra in slit position S2 are extracted from the 
circumnuclear regions located to the North and North- West 
of the nucleus, and therefore any contribution from the un- 
derlying galaxy bulge is difficult to assess. Slit position SI 
crosses the galactic nucleus. This can be used to estimate the 
underlying bulge contribution. For the blue spectra, it turns 
out to be almost negligible amounting to, at most, 10 per 
cent at the H/3 line. For the red spectra, the bulge contribu- 
tion is more important. From Gaussian fits to the A 8620 A 
continuum profile of SI we find it to be about 25 per cent 
for R4, the weakest region, and the one closest to the nu- 
cleus. On the other hand, the analysis of the broad near-IR 
HST-NICMOS images shown in Fig. 1 shows less contrast 
between the emission from the regions and the underlying 
bulge which is very close to the image background emission. 
Its contribution is about 25 per cent for the weak regions in 
the central zone of NGC 3310, that is in very good agreement 
with the cluster identification made by Elmegreen et al. 
(2002) using the equivalent J and K-band HST-NICMOS 
images and the ground base data from KPNO (J and Ri- 
bands) . 

Fig. 5 shows the spectra of the observed circumnuclear 
regions split into two panels corresponding to the blue and 
the red spectral ranges. The spectrum of the nucleus of 
NGC 3310 is shown in Fig. 6. The blue spectra show the 
Balmer H/3 recombination line and the collisionally excited 



[Olll] lines at AA 4959, 5007 A. Generally, these forbidden 
lines are very weak (see for example Figs. 5 of Papers I 
and II), and, in some cases, only the strongest A 5007 A is 
detected (right-hand panels of these Figs.). However, in the 
case of NGC 3310 they are very strong, due to the low abun- 
dance of the CNSFRs in this galaxy, with values between 
0.2-0.4 Z (Pastoriza et al. 1993). These low values of the 
abundances can be explained by the probably unusual inter- 
action history of the galaxy (Elmegreen et al. 2002; Rregel 
& Sancisi 2001; Mulder & van Driel 1996; Smith et al. 1996; 
Balick & Heckman 1981), fuelling the ring with accreted 
neutral gas, as modeled by Athanassoula (1992) and Piner 
et al. (1995). The red spectra show the stellar CaT lines in 
absorption. In some cases these lines are contaminated by 
Paschen emission which occurs at wavelengths very close 
to those of the CaT lines. Other emission features, such 
as OiA8446, [Cln]A8579, Pa 14 and [Fen] A 8617 are also 
present. In Fig. 5 we can easily appreciate, for example in 
R4+R5, the Paschen series from Pal 3 to Pa22, as well as 
the previously mentioned lines of O, CI and Fe. In all cases, 
a single Gaussian fit to the emission lines was performed 
and the lines were subsequently subtracted (see also Ostlin 
et al. 2004; Cumming et al. 2008) after checking that the 
theoretically expected ratio between the Paschen lines was 
satisfied. The observed red spectra are plotted with a dashed 
line. The solid line shows the subtracted spectra. 

Fig. 7 shows the spectra of the almost pure emission 
knots labelled X and Y, and those of the Jumbo region. In 
all cases the blue range of the spectrum presents very intense 
emission lines. The red spectral range presents a very weak 
and noisy continuum. In the case of region X only noise is 
detected therefore no spectrum is shown in Fig. 7. In the 
other two regions we can see a set of emission lines. The 
Jumbo region presents many strong lines, even He and Ni in 
emission. Due to the low signal-to-noise ratio of the contin- 
uum and the presence of the strong emission lines in these 
regions of NGC 3310, we could not obtain stellar spectra 
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Figure 3. Spatial profiles of H/3, [Olll] A5007A and the far red (upper, middle and lower panels respectively) for each slit. For the 
emission lines, the profiles correspond to line+continuum (dashed line), continuum (dashed-dotted line) and the difference between them 
(solid line), representing the pure emission from Yi(3 and [Olll] respectively. For the far red profiles, the solid lines represent the continuum. 
Pixel number increases to the North. Horizontal small lines show the location of the CNSFRs and nuclear apertures. 
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Figure 4. Enlargement of the spatial profiles presented in Fig. 3, except in the far red range of S2. 
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Figure 5. Blue (left) and red (right) rest frame normalized spectra of the observed CNSFRs. For the red range, the dashed line shows 
the observed spectrum; the solid line represents the spectrum after subtracting the emission lines (see text). 
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Figure 5. (cont) Blue (left) and red (right) rest frame normalized spectra of the observed CNSFRs. For the red range, the dashed line 
shows the observed spectrum; the solid line represents the spectrum after subtracting the emission lines (see text). 



with enough signal in the CaT absorption feature to allow 
an accurate measurement of velocity dispersions. 



3.1 Kinematics of stars and ionized gas 

A detailed description of the methods and techniques used to 
derive the values of radial velocities and velocity dispersions 
as well as sizes, masses and emission line fluxes has been 
given in Paper I. Therefore only a brief summary is given 
below. 



Stellar analysis 

Stellar radial velocities and velocity dispersions were ob- 
tained from the CaT absorption lines using the cross- 
correlation technique, described in detail by Tonry & Davis 
(1979). This method requires the comparison with a stellar 
template that represents the stellar population that best re- 
produces the absorption features. This has been built from 
a set of 11 late- type giant and supergiant stars with strong 
CaT absorption lines. We have followed the work by Nelson 
& Whittle (1995) with the variation introduced by Palacios 
et al. (1997) of using the individual stellar templates in- 
stead of an average. This procedure will allow us to correct 
for the known possible mismatches between template stars 
and the region composite spectrum. The implementation of 
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the method in the external package of iraf XCSAO (Kurtz 
& Mink 1998) has been used. 

To determine the line-of-sight stellar velocity and veloc- 
ity dispersion along each slit, extractions were made every 
two pixels for slit position SI and every three pixels for slit 
position S2, with one pixel overlap between consecutive ex- 
tractions in this latter case. In this way the S/N ratio and 
the spatial resolution were optimized. Besides, the stellar 
velocity dispersion was estimated at the position of each 
CNSFR and the nucleus using an aperture of five pixels in 
all cases, which corresponds to 1.0 x 1.8arcsec 2 . The veloc- 
ity dispersion (cr) of the stars (cr*) is taken as the average of 
the cr values found for each stellar template, and its error is 
taken as the dispersion of the individual values of a and the 
rms of the residuals of the wavelength fit. These values are 
listed in column 3 of Table 2 along with their corresponding 
errors. Stellar velocity dispersions of X, Y and the Jumbo 
region could not be estimated due to the low signal-to-noise 
ratio of the continua and the CaT absorption features. The 
same is true for region R1+R2, where the red continuum 
and the CaT features after subtracting the emission lines 
have a low signal-to- noise ratio, although an estimate of the 
stellar velocity dispersion could be given in this case. 

The radial velocities have been determined directly from 
the position of the main peak of the cross-correlation of each 
galaxy spectrum with each template in the rest frame. The 
average of these values is the final adopted radial velocity. 

Ionized gas analysis 

The velocity dispersion of the ionized gas was estimated 
for each observed CNSFR and for the galaxy nucleus from 
Gaussian fits to the K/3 and [Olll] A 5007 A emission lines us- 
ing five pixel apertures, corresponding to 1.0 x 1.9arcsec 2 . 
For a single Gaussian fit, the position and width of a given 
emission line is taken as the average of the fitted Gaus- 
sians to the whole line using three different suitable continua 
(Jimenez-Benito et al. 2000), and their errors are given by 
the dispersion of these measurements taking into account the 
rms of the wavelength calibration. In all the studied regions, 
however, the best fit for the emission lines are obtained with 
two different components having different radial velocities of 
up to 25 km s _1 . The radial velocities found for the narrow 



and broad components of both H/3 and [Olll] , are the same 
within the errors. An example of the two-Gaussian fit for 
R4+R5 in S2 is shown in Fig. 8. 

For each CNSFR, the gas velocity dispersion for the H[3 
and [Olll] A 5007 A lines derived using single and double line 
Gaussian fits, and their corresponding errors are listed in 
Table 2. Columns 4 and 5, labelled 'One component', give 
the results for the single Gaussian fit. Columns 6 and 7, 
and 8 and 9, labelled 'Two components - Narrow' and 'Two 
components - Broad', respectively, list the results for the two 
component fits. The last column of the table, labelled Av n b, 
gives the velocity difference between the narrow and broad 
components. This is calculated as the average of the H/3 and 
[Olll] fit differences. Taking into account the errors in the 
two component fits, the errors in these velocity differences 
vary from 5 to lOkrns -1 . 

We have also determined the distribution along each slit 
position of the radial velocities and the velocity dispersions 
of the ionized gas using the same procedure as for the stars, 
that is using spectra extracted every two pixels for SI and 
every three pixels, superposing one pixel for consecutive ex- 
tractions, for S2. These spectra, however, do not have the 
required S/N ratio to allow an acceptable two-component fit, 
therefore a single- Gaussian component has been used. The 
goodness of this procedure is discussed in Section 4 below. 

3.2 Emission line ratios 

We have used two different ways to integrate the intensity 
of a given line: (1) if an adequate fit was attained by a single 
Gaussian, the emission line intensities were measured using 
the splot task in IRAF. For the H/3 emission lines a conspic- 
uous underlying stellar population is inferred from the pres- 
ence of absorption features that depress the lines (e.g. see 
discussion in Diaz 1988). Examples of this effect can be ap- 
preciated in Fig. 9. We have defined a pseudo-continuum at 
the base of the line to measure the line intensities and min- 
imize the errors introduced by the underlying population 
(for details see Hagele et al. 2006). (2) When the optimal fit 
was obtained by two Gaussians the individual intensities of 
the narrow and broad components are estimated from the 
fitting parameters (I = 1.0645 A x FWHM = a/2tt A cr; where 
I is the Gaussian intensity, A is the amplitude of the Gaus- 
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Figure 7. Idem as Fig. 5 for regions J, X (only blue) and Y. The dashed line shows the observed spectrum; the solid line represents the 
spectrum after subtracting the emission lines (see text). 



sian, FWHM is the full width at half- maximum and a is the 
dispersion of the Gaussian). A pseudo-continuum for the 
11(3 emission line was also defined in these cases. The sta- 
tistical errors associated with the observed emission fluxes 
have been calculated with the expression ai = cr c N 1//2 [l + 
EW/(NA)] 1/2 , where ai is the error in the observed line 
flux, a c represents the standard deviation in a box near the 
measured emission line and stands for the error in the con- 
tinuum placement, N is the number of pixels used in the 
measurement of the line intensity, EW is the line equiva- 
lent width and A is the wavelength dispersion in A pixel -1 
(Gonzalez-Delgado et al. 1994). For the H/3 emission line 
we have doubled the derived error, ai, in order to take into 



account the uncertainties introduced by the presence of the 
underlying stellar population (Hagele et al. 2006). 

The logarithmic ratio between the emission line intensi- 
ties of [Oiii] A 5007 A and H/3 and their corresponding errors 
are presented in Table 3. We have also listed the logarithmic 
ratio between the emission line fluxes of [Nil] A 6584 A and 
Ha together with their corresponding errors. These values 
listed for R4, R4+R5, S6, J and the nucleus are from Pastor- 
iza et al. (1993; S6 seems to be their region L). For the rest 
of the regions we used an extrapolation of the results given 
by them from the spatial profiles of these emission lines. 
They reported a constant value of 0.2 for the [Nil] / Ha ra- 
tio for the Hn regions in their slit positions 2 and 3, while 
this value changes from 0.2 to 0.5 along position 1 over the 
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Table 2. Velocity dispersions. 



1 component 2 components 

narrow broad 



Region 


Slit 


(7* 


0-gas(W) 


o- g as([Oin]) 


0-gas(tt{3) 


^ os ([Oiii]) 


0-gas(W) 


°gas{[Om]) 


Av n6 


R1+R2 


S2 


80: 


33±4 


31±4 


24±5 


22±5 


54±7 


50±4 


-10 


R4 


SI 


36±3 


34±3 


32±3 


28±4 


26±3 


55±5 


52±4 


10 


R4+R5 


S2 


38±3 


27±4 


22±3 


22±4 


18±3 


46±5 


40±2 


15 


R6 


S2 


35±5 


30±3 


28±3 


23±3 


21±3 


54±7 


56±3 





S6 


S2 


31±4 


27±3 


26±3 


20±3 


19±3 


47±4 


47±3 


10 


R7 


S2 


44±5 


21±5 


17±4 


18±4 


14±3 


41±4 


36±3 





RIO 


SI 


39±3 


38±3 


40±3 


26±2 


26±3 


54±2 


59±4 


-20 


N 


SI 


73±3 


55±3 


66±3 


35±3 


35±3 


73±3 


83±4 


20 


J 


SI 




34±2 


30±2 


25±2 


22±2 


61±3 


57±3 


-25 


X 


S2 




22±4 


18±3 


18±4 


14±3 


40±5 


30±3 


-5 


Y 


SI 




28±3 


28±3 


22±4 


25±3 


44±4 


61±5 


10 



velocity dispersions in kms 




Figure 8. Sections of the normalized spectrum of R4+R5 (solid line). The left panel shows from 4868 to 4883 A containing ¥Lf3 and 
the right panel shows from 5014 to 5029 A containing the [Om] A 5007 A emission line. For both we have superposed the fits from the 
NGAUSSFIT task in IRAF; the dashed-dotted line is the broad component, the dotted line is the narrow component and the dashed line is 
the sum of both. 



nucleus of the galaxy, being 0.23 and 0.27 for regions B and 
L, respectively. We adopt for R1+R2, R6, R7, R10, X and 
Y, a value of 0.2 without assigning any error to it. 



4 DYNAMICAL MASS DERIVATION 

The mass of a virialized stellar system is given by two pa- 
rameters: its velocity dispersion and its size. 

In order to determine the sizes of the stellar clusters 
within our observed CNSFRs, we have used the retrieved 
wide V HST image which provides a spatial resolution of 
0.045 " per pixel. At the distance of NGC 3310, this resolu- 
tion corresponds to 3.3pc/pixel. Fig. 10 shows enlargements 
around the different regions studied with intensity contours 
overlapped. 

We find, as expected, that the eight CNSFRs stud- 
ied here are formed by a large number of individual star- 
forming clusters. When we analyze the F606W-PC1 image 
of NGC 3310, we find another region very close to R6, not 
classified by Diaz et al. (2000), probably due to the lower 
spatial resolution of the data used in that work, and la- 
belled S6 by us, which seems to be coincident with knot L 
of Pastoriza et al. (1993). For all the regions of this galaxy 



we find a principal knot and several secondary knots with 
lower peak intensities, except for Rl and R2, which present 
only one and two secondary knots, respectively. The same 
criterion as in Papers I and II has been applied to name the 
regions. For the other regions of this galaxy, R4, R5, R6, 
S6, R7 and R10, we find 31, 8, 14, 9, 29 and 15 individ- 
ual clusters, respectively. In all cases the knots have been 
found with a detection level of 10 a above the background. 
All these knots are within the radius of the regions denned 
by Diaz et al. (2000), except for S6. We have to remark that 
our search for knots has not been exhaustive since that is 
not the aim of this work. 

We have fitted circular regions to the intensity contours 
corresponding to the half light brightness distribution of 
each single structure (see Fig. 10), following the procedure 
given in Meurer et al. (1995), assuming that the regions have 
a circularly symmetric Gaussian profile. The radii of the sin- 
gle knots vary between 2.2 and 6.2 pc. Table 4 gives, for each 
identified knot, the position, as given by the astrometric cal- 
ibration of the HST image; the radius of the circular region 
denned as described above together with its error; and the 
peak intensity in counts, as measured from the WFPC2 im- 
age. For the nucleus of NGC 3310 we measure a radius of 
14.2 pc. 
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Figure 9. Enlargement of the blue rest frame normalized spectra of R4 (left) and R6 (right). 



Table 3. Line ratios. 



One component 



Two components 



Region 


Slit 


log([Oin]5007/H/?) 


Narrow 
log([Oin]5007/H/?) 


Broad 
log([Oin]5007/H/?) 


log([Nn]6584/Ha) a 


R1+R2 


S2 


0.42±0.01 


0.39±0.01 


0.45±0.03 


-0.70: 


R4 


SI 


0.28±0.01 


0.28±0.01 


0.28±0.04 


-0.69±0.01 


R4+R5 


S2 


0.42±0.01 


0.39±0.01 


0.44±0.03 


-0.69±0.01 


R6 


S2 


0.28±0.01 


0.30±0.01 


0.28±0.06 


-0.70: 


S6 


S2 


0.21±0.01 


0.21±0.01 


0.22±0.08 


-0.56±0.01 


R7 


S2 


0.23±0.01 


0.21±0.01 


0.32±0.09 


-0.70: 


RIO 


SI 


0.19±0.01 


0.20±0.04 


0.20±0.08 


-0.70: 


N 


S2 


0.04±0.02 


-0.13±0.06 


0.11±0.07 


-0.30±0.01 


J 


SI 


0.45±0.01 


0.45±0.01 


0.45±0.01 


-0.80±0.01 


X 


S2 


0.60±0.01 


0.55±0.01 


0.67±0.04 


-0.70: 


Y 


SI 


0.28±0.01 


0.43±0.01 


0.00±0.04 


-0.70: 


a From Pastoriza et al. (1993). 



Upper limits to the dynamical masses (M*) inside the 
half light radius (R) for each observed knot have been 
estimated under the following assumptions: (i) the sys- 
tems are spherically symmetric; (ii) they are gravitationally 
bound; and (iii) they have isotropic velocity distributions 
[a 2 (total) = 3 a 2 ]. The general expression for the virial mass 
of a cluster is r\ a 2 R/G, where R is the effective gravita- 
tional radius and 77 is a dimensionless number that takes 
into account departures from isotropy in the velocity dis- 
tribution and the spatial mass distribution, binary fraction, 
mean surface density, etc. (Boily et al. 2005; Fleck et al. 
2006). Following Ho & Filippenko (1996a,b), and for consis- 
tence with Paper I, Paper II and Hagele (2008), we obtain 
the dynamical masses inside the half-light radius using 77 = 3 
and adopting the half-light radius as a reasonable approx- 
imation of the effective radius. Other authors (e.g. Spitzer 
1987; Smith & Gallagher 2001; Moll et al. 2008) assumed 
that the 77 value is about 9.75 to obtain the total mass. 
On the absence of any knowledge about the tidal radius of 
the clusters, we adopted this conservative approach. On the 
derived masses, the different adopted 77 values act as multi- 
plicative factors. 

It must be noted that while we can measure the size of 
each knot, we do not have direct access to the stellar velocity 
dispersion of each individual cluster, since our spectroscopic 
measurements encompass a wider area (1.0 x 1.9 arcsec 2 
which corresponds approximately to 73 x 131 pc 2 at the 



adopted distance for NGC 3310) that includes the whole 
CNSFRs to which each group of knots belongs. The use of 
these wider size scale velocity dispersion measurements to 
estimate the mass of each knot, leads us to overestimate the 
mass of the individual clusters, and hence of each CNSFR. 
As we can not be sure that we are actually measuring their 
velocity dispersion, we prefer to say that our measurements 
of a*, and hence the dynamical masses, constitute upper 
limits. Although we are well aware of the difficulties, still we 
are confident that these upper limits are valid and important 
for comparison with the gas kinematical measurements. 



The estimated dynamical masses for each knot and their 
corresponding errors are listed in Table 5. For the regions 
that have been observed in more than one slit position, we 
list the derived values using the two separate stellar veloc- 
ity dispersions. The dynamical masses in the rows labelled 
"sum" have been found by adding the individual masses 
in a given CNSFR, as well as the galaxy nucleus, N. The 
fractional errors of the dynamical masses of the individual 
knots and of the CNSFRs are listed in column 4. As ex- 
plained above, since the stellar velocity dispersion for the 
R1+R2 region of NGC 3310 has a large error, we do not list 
the errors of the dynamical masses of the individual knots 
or of the whole region R12sum. 
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Table 4. Positions, radii and peak intensities derived from the HST-F606W image. 



Region 


position 




R 


I 


Region 


position 


R 


I 




^72000.0 


<5j2000.0 




(P<0 


(counts) 




^72000.0 


<5j2000.0 


(PC) 


(counts) 


Rl 


10 n 38 m 47f070 


+53° 30' 14'. 


'46 


4.6±0.3 


322 


S6 


10 n 38 m 45!748 


+53°30'18 / /42 


4.2±0.2 


544 


Rla 


10 h 38 m 47?032 


+53° 30' 14'. 


'74 


2.6±0.2 


179 


S6a 


10 h 38 m 45!768 


+53 o 30'18 / /00 


5.5±0.2 


385 


R2 


10 h 38 m 47f026 


+53°30 / 15 / . 


'55 


4.7±0.2 


1271 


S6b 


10 h 38 m 45!700 


+53°30'18 / /26 


2.6±0.3 


237 


R2a 


10 h 38 m 46?997 


+53°30 / 15 / . 


'52 


5.5±0.5 


515 


S6c 


10 h 38 m 45!711 


+53 o 30'18 / /14 


3.4±0.3 


204 


R2b 


10 h 38 m 47f008 


+53°30'15'. 


'84 


3.3±0.2 


287 


S6d 


10 h 38 m 45!813 


+53°30'17 / /86 


2.7±0.3 


180 


R4 


10 h 38 m 46!109 


+53°30 / 16 / . 


'04 


2.9±0.2 


2409 


S6e 


10 h 38 m 45!770 


+53 o 30'18 / /28 


3.0±0.2 


174 


R4a 


10 h 38 m 46!198 


+53°30 / 17 / . 


'01 


3.1±0.2 


1942 


S6f 


10 h 38 m 45!842 


+53 o 30'18 / /00 


3.6±0.2 


173 


R4b 


10 h 38 m 46!242 


+53°30 / 16 / . 


'57 


2.5±0.2 


1038 


S6g 


10 h 38 m 45!801 


+53 o 30'18 / /02 


2.9±0.3 


172 


R4c 


10 h 38 m 46!311 


+53°30 / 15 / . 


'66 


2.4±0.2 


584 


S6h 


10 h 38 m 45!738 


+53 o 30'17 / /98 


2.8±0.3 


166 


R4d 


10 h 38 m 46!115 


+53°30 / 16 / . 


'25 


4.1±0.2 


494 


R7 


10 h 38 m 45!270 


+53 o 30'18 / /18 


3.1±0.2 


639 


R4e 


10 h 38 m 46!319 


+53°30 / 16 / . 


'67 


4.2±0.3 


489 


R7a 


10 h 38 m 45!520 


+53 o 30'18 / /24 


5.3±0.2 


599 


R4f 


10 h 38 m 46!188 


+53°30 / 16 / . 


'59 


4.2±0.2 


403 


R7b 


10 h 38 m 45!537 


+53 o 30'18 / /16 


3.5±0.3 


449 


R4g 


10 h 38 m 46!264 


+53°30 / 17 / . 


'81 


2.9±0.2 


391 


R7c 


10 h 38 m 45!390 


+53°30 / 19 / /14 


4.8±0.2 


390 


R4h 


10 h 38 m 46!210 


+53°30 / 16 / . 


'47 


2.8±0.4 


336 


R7d 


10 h 38 m 45!633 


+53 o 30'17 / /57 


3.0±0.2 


275 


R4i 


10 h 38 m 46!353 


+53° 30' 14'. 


'68 


4.2±0.2 


304 


R7e 


10 h 38 m 45!527 


+53 o 30'18 / /86 


2.4±0.3 


273 


R4j 


10 h 38 m 46!274 


+53°30'16 / . 


'36 


3.0±0.3 


284 


R7f 


10 h 38 m 45!466 


+53 o 30'18 / /03 


3.9±0.2 


264 


R4k 


10 h 38 m 46!107 


+53°30 / 15 / . 


'56 


2.7±0.2 


271 


R7g 


10 h 38 m 45!607 


+53 o 30'17 / /05 


4.5±0.2 


231 


R41 


10 h 38 m 46!258 


+53°30 / 15 / . 


'30 


3.4±0.5 


268 


R7h 


10 h 38 m 45!613 


+53 o 30'19 / /21 


3.9±0.2 


209 


R4m 


10 h 38 m 46!513 


+53°30'15'. 


'02 


2.8±0.3 


259 


R7i 


10 h 38 m 45!394 


+53 o 30'17 / /51 


3.0±0.2 


201 


R4n 


10 h 38 m 46!316 


+53°30 / 17 / . 


'66 


3.4±0.3 


250 


R7j 


10 h 38 m 45!638 


+53 o 30'17 / /45 


4.0±0.5 


200 


R4o 


10 h 38 m 46!236 


+53°30 / 15 / . 


'69 


3.1±0.2 


241 


R7k 


10 h 38 m 45!313 


+53 o 30'17 / /00 


3.4±0.3 


198 


R4p 


10 h 38 m 46!327 


+53° 30' 14'. 


'82 


2.4±0.2 


241 


R71 


10 h 38 m 45!371 


+53 o 30'17 / /53 


4.1±0.3 


191 


R4q 


10 h 38 m 46!352 


+53° 30' 14'. 


'53 


3.1±0.4 


227 


R7m 


10 h 38 m 45!603 


+53°30 / 18 / /39 


3.0±0.3 


190 


R4r 


10 h 38 m 46!200 


+53°30 / 16 / . 


'77 


3.1±0.3 


216 


R7n 


10 h 38 m 45!363 


+53 o 30'17 / /10 


3.4±0.3 


180 


R4s 


10 h 38 m 46!279 


+53°30 / 16 / . 


'72 


2.7±0.3 


212 


R7o 


10 h 38 m 45!549 


+53 o 30'17 / /87 


3.9±0.5 


179 


R4t 


10 h 38 m 46!123 


+53°30'16'. 


'40 


2.9±0.4 


206 


R7p 


10 h 38 m 45!461 


+53°30'18 / /84 


3.2±0.2 


177 


R4u 


10 h 38 m 46?350 


+53° 30' 14'. 


'93 


3.8±0.4 


206 


R7q 


10 h 38 m 45!494 


+53 o 30'16 / /98 


3.8±0.3 


170 


R4v 


10 h 38 m 46!347 


+53°30 / 15 / . 


'14 


4.0±0.4 


205 


R7r 


10 h 38 m 45!533 


+53 o 30'16 / /82 


3.9±0.3 


168 


R4w 


10 h 38 m 46!291 


+53°30 / 16 / . 


'18 


2.9±0.4 


204 


R7s 


10 h 38 m 45!576 


+53 o 30'18 / /10 


4.0±0.2 


165 


R4x 


10 h 38 m 46!308 


+53°30'15'. 


'19 


3.1±0.3 


201 


R7t 


10 h 38 m 45!422 


+53 o 30'18 / /27 


2.4±0.3 


162 


R4y 


10 h 38 m 46!313 


+53°30 / 16 / . 


'20 


2.9±0.4 


194 


R7u 


10 h 38 m 45!306 


+53 o 30'17 / /08 


3.7±0.4 


160 


R4z 


10 h 38 m 46!366 


+53° 30' 14'. 


'81 


3.5±0.3 


191 


R7v 


10 h 38 m 45!364 


+53 o 30'18 / /99 


3.7±0.3 


157 


R4a 


10 h 38 m 46!282 


+53°30 / 16 / . 


'48 


3.3±0.4 


187 


R7w 


10 h 38 m 45!423 


+53 o 30'19 / /19 


3.9±0.2 


154 


R4(3 


10 h 38 m 46!296 


+53°30'15'. 


'04 


3.7±0.3 


183 


R7x 


10 h 38 m 45!447 


+53°30'18 / /24 


3.5±0.2 


152 


R4 7 


10 h 38 m 46!140 


+53°30 / 16 / . 


'19 


3.7±0.3 


166 


R7y 


10 h 38 m 45!550 


+53 o 30'16 / /96 


3.9±0.3 


147 


R4S 


10 h 38 m 46!136 


+53°30'15'. 


'94 


2.6±0.3 


166 


R7z 


10 h 38 m 45!476 


+53 o 30'18 / /34 


3.3±0.4 


145 


R5 


10 h 38 m 46!lll 


+53°30 / 17 / . 


'03 


3.4±0.3 


606 


R7a 


10 h 38 m 45!543 


+53 o 30'18 / /70 


4.7±0.6 


143 


R5a 


10 h 38 m 46!119 


+53°30'17'. 


'42 


3.1±0.4 


370 


R7(3 


10 h 38 m 45!545 


+53°30'18 / /99 


3.4±0.4 


138 


R5b 


10 h 38 m 46!094 


+53°30 / 17 / . 


'26 


3.6±0.4 


313 


R10 


10 h 38 m 45!309 


+53 o 30'08 / /14 


3.1±0.2 


354 


R5c 


10 h 38 m 46!104 


+53°30 / 17 / . 


'79 


6.2±0.6 


262 


RlOa 


10 h 38 m 45!335 


+53 o 30'07 / /59 


5.1±0.5 


303 


R5d 


10 h 38 m 46?153 


+53°30 / 17 / . 


'08 


3.6±0.3 


244 


RlOb 


10 h 38 m 45!420 


+53°30'08 / /06 


3.8±0.2 


270 


R5e 


10 h 38 m 46!071 


+53°30'16'. 


'97 


3.7±0.2 


222 


RIOc 


10 h 38 m 45!343 


+53 o 30'07 / /73 


3.5±0.3 


260 


R5f 


10 h 38 m 46?054 


+53°30 / 17 / . 


'48 


3.0±0.2 


206 


RIOd 


10 h 38 m 45!338 


+53 o 30'07 / /84 


4.1±0.4 


251 


R5g 


10 h 38 m 46!083 


+53°30 / 16 / . 


'71 


3.2±0.3 


193 


RlOe 


10 h 38 m 45!307 


+53 o 30'08 / /28 


2.8±0.3 


235 


R6 


10 h 38 m 45?935 


+53°30 / 18 / . 


'20 


3.4±0.2 


378 


RIOf 


10 h 38 m 45!298 


+53 o 30'08 / /41 


2.8±0.4 


222 


R6a 


10 h 38 m 46!009 


+53°30'17'. 


'83 


2.2±0.3 


342 


RlOg 


10 h 38 m 45!393 


+53 o 30'07 / /84 


4.4±0.6 


200 


R6b 


10 h 38 m 45!879 


+53°30 / 17 / . 


'69 


3.9±0.2 


260 


RIOh 


10 h 38 m 45!254 


+53 o 30'07 / /38 


4.1±0.3 


199 


R6c 


10 h 38 m 45!972 


+53°30 / 17 / . 


'51 


2.7±0.2 


258 


RIOi 


10 h 38 m 45!395 


+53 o 30'07 / /70 


4.7±0.5 


182 


R6d 


10 h 38 m 45!907 


+53°30 / 17 / . 


'88 


4.2±0.3 


242 


RIOj 


10 h 38 m 45!371 


+53°30 / 08'/15 


4.0±0.5 


177 


R6e 


10 h 38 m 46f005 


+53°30'17'. 


'46 


3.9±0.4 


221 


RIOk 


10 h 38 m 45!421 


+53°30'08 / /27 


3.8±0.3 


177 


R6f 


10 h 38 m 45?952 


+53°30 / 18 / . 


'44 


3.7±0.2 


215 


R101 


10 h 38 m 45!408 


+53 o 30'07 / /80 


3.9±0.3 


173 


R6g 


10 h 38 m 45!907 


+53°30 / 16 / . 


'99 


3.4±0.4 


170 


RIOm 


10 h 38 m 45!374 


+53 o 30'07 / /86 


3.1±0.3 


169 


R6h 


10 h 38 m 45!948 


+53°30'17 / . 


'94 


3.1±0.3 


168 


RIOn 


10 h 38 m 45!424 


+53 o 30'07 / /90 


3.1±0.4 


168 


R6i 


10 h 38 m 45!965 


+53°30 / 17 / . 


'75 


2.5±0.2 


168 


N 


10 h 38 m 45!893 


+53 o 30'll / /47 


14.2±1.0 


3735 


R6j 


10 h 38 m 46!007 


+53°30 / 16 / . 


'94 


4.0±0.5 


165 


Na 


10 h 38 m 45!970 


+53 o 30'10 / /85 


5.8±0.6 


370 


R6k 


10 h 38 m 45!988 


+53°30 / 17 / . 


'28 


3.1±0.2 


160 












R61 


10 h 38 m 46!001 


+53°30 / 17 / . 


'19 


3.7±0.4 


156 












R6m 


10 h 38 m 45!967 


+53°30 / 17 / . 


'12 


4.2±0.5 


146 
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5 IONIZING STAR CLUSTER PROPERTIES 

For each of the CNSFR the total number of ionizing pho- 
tons was derived from the total observed Ha luminosities 
given by Diaz et al. (2000) and Pastoriza et al. (1993), cor- 
recting for the different assumed distance. Diaz et al. (2000) 
and Pastoriza et al. (1993) give the already extinction cor- 
rected luminosities. Diaz et al. (2000) estimated a diameter 
of 2arcsec for regions R2, R4, R6, R7 and R10, 1.4arcsec 
for Rl and R5, and 3.4 for the Jumbo region. For this latter 



region we give the values derived using the different quan- 
tities for the luminosities given by (Diaz et al. 2000, region 
R19) and (Pastoriza et al. 1993, region A). For R1+R2 and 
R4+R5 (S2) we added their Ha luminosities. No values are 
found in the literature for the Ha luminosity of regions X 
and Y. Our derived values of Q(Ho) constitute lower limits 
since we have not taken into account the fraction of photons 
that may have been absorbed by dust or may have escaped 
the region. 

Once the number of Lyman continuum photons have 
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Table 5. Upper limits to the dynamical masses. 



Region 


Slit 


M* 




error(%) 


Region 


Slit 


M* 


error(%) 


Region 


Slit 


M* 


error(%) 


Rl 


S2 


202 







R4n 


S2 


35±7 


19 


R7 


S2 


42±9 


22 


Rla 


S2 


115 







R4o 


S2 


32±6 


18 


R7a 


S2 


71±15 


21 


R2 


S2 


208 







R4p 


S2 


25±5 


19 


R7b 


S2 


46±11 


23 


R2a 


S2 


241 







R4q 


S2 


32±7 


21 


R7c 


S2 


64±14 


21 


R2b 


S2 


144 







R4r 


S2 


31±6 


19 


R7d 


S2 


40±9 


22 


R12sum 


S2 


911 







R4s 


S2 


27±6 


20 


R7e 


S2 


32±8 


24 


R4 


SI 


26±4 


17 


R4t 


S2 


30±6 


22 


R7f 


S2 


52±11 


22 


R4a 


SI 


27±5 


17 


R4u 


S2 


39±8 


20 


R7g 


S2 


60±13 


22 


R4b 


SI 


22±4 


18 


R4v 


S2 


41±8 


19 


R7h 


S2 


51±11 


22 


R4c 


SI 


21±4 


18 


R4w 


S2 


30±6 


22 


R7i 


S2 


40±9 


22 


R4d 


SI 


36±6 


16 


R4x 


S2 


32±6 


19 


R7j 


S2 


53±13 


24 


R4e 


SI 


38±6 


17 


R4y 


S2 


30±6 


22 


R7k 


S2 


44±10 


23 


R4f 


SI 


38±6 


16 


R4z 


S2 


36±7 


19 


R71 


S2 


54±12 


22 


R4g 


SI 


26±4 


17 


R4a 


S2 


33±7 


21 


R7m 


S2 


40±9 


23 


R4h 


SI 


25±5 


21 


R4/3 


S2 


38±7 


19 


R7n 


S2 


44±10 


23 


R4i 


SI 


37±6 


16 


R4 7 


S2 


38±7 


19 


R7o 


S2 


52±13 


25 


R4j 


SI 


27±5 


19 


R45 


S2 


27±5 


20 


R7p 


S2 


43±9 


22 


R4k 


SI 


24±4 


17 


R5 


S2 


35±7 


19 


R7q 


S2 


50±11 


22 


R41 


SI 


30±7 


21 


R5a 


S2 


32±7 


21 


R7r 


S2 


52±12 


22 


R4m 


SI 


25±5 


19 


R5b 


S2 


36±7 


20 


R7s 


S2 


53±12 


22 


R4n 


SI 


30±5 


18 


R5c 


S2 


63±12 


19 


R7t 


S2 


32±8 


24 


R4o 


SI 


28±5 


17 


R5d 


S2 


36±7 


19 


R7u 


S2 


48±11 


24 


R4p 


SI 


21±4 


18 


R5e 


S2 


37±7 


18 


R7v 


S2 


48±11 


23 


R4q 


SI 


28±6 


20 


R5f 


S2 


30±5 


18 


R7w 


S2 


52±11 


22 


R4r 


SI 


27±5 


18 


R5g 


S2 


33±6 


19 


R7x 


S2 


46±10 


22 


R4s 


SI 


24±5 


19 


R45sum 


S2 


1317±41 


3 


R7y 


S2 


52±12 


22 


R4t 


SI 


26±5 


21 


R6 


S2 


28±8 


30 


R7z 


S2 


44±11 


24 


R4u 


SI 


34±6 


19 


R6a 


S2 


18±6 


32 


R7a 


S2 


63±15 


25 


R4v 


SI 


36±7 


19 


R6b 


S2 


32±9 


29 


R7^ 


S2 


44dzll 


24 


R4w 


SI 


26±5 


21 


R6c 


S2 


22±7 


30 


R7sum 


S2 


1413±60 


4 


R4x 


SI 


28±5 


18 


R6d 


S2 


35±10 


30 


RIO 


SI 


33±5 


17 


R4y 


SI 


26±5 


21 


R6e 


S2 


32±10 


31 


RlOa 


SI 


53±10 


18 


R4z 


SI 


31±6 


18 


R6f 


S2 


30±9 


29 


RlOb 


SI 


40±6 


16 


R4a 


SI 


29±6 


20 


R6g 


S2 


28±9 


31 


RIOc 


SI 


37±6 


18 


R4/3 


SI 


33±6 


18 


R6h 


S2 


26±8 


30 


RIOd 


SI 


43±8 


18 


R4 7 


SI 


33±6 


18 


R6i 


S2 


21±6 


30 


RlOe 


SI 


29±5 


19 


R45 


SI 


23±5 


19 


R6j 


S2 


33±10 


31 


RIOf 


SI 


29±6 


21 


R4sum 


SI 


886±30 


3 


R6k 


S2 


25±8 


30 


RlOg 


SI 


46±9 


21 


R4 


S2 


30±5 


18 


R61 


S2 


31±10 


31 


RIOh 


SI 


43±7 


17 


R4a 


S2 


31±6 


18 


R6m 


S2 


35dzll 


31 


RIOi 


SI 


50±9 


19 


R4b 


S2 


25±5 


19 


R6sum 


S2 


397±33 


8 


RIOj 


SI 


42±8 


19 


R4c 


S2 


25±5 


19 


S6 


S2 


28±6 


23 


RIOk 


SI 


40±7 


17 


R4d 


S2 


42±7 


17 


S6a 


S2 


37±8 


23 


R101 


SI 


40±7 


17 


R4e 


S2 


43±8 


18 


S6b 


S2 


18dz5 


25 


RIOm 


SI 


32±6 


18 


R4f 


S2 


43±7 


17 


S6c 


S2 


23±6 


24 


RIOn 


SI 


33±7 


20 


R4g 


S2 


30±5 


18 


S6d 


S2 


18±5 


25 


RIOsum 


SI 


588±28 


5 


R4h 


S2 


28±6 


22 


S6e 


S2 


20±5 


23 


N 


SI 


526±57 


11 


R4i 


S2 


42±7 


17 


S6f 


S2 


24±6 


23 


Na 


SI 


216±28 


13 


R4j 


S2 


30±6 


19 


S6g 


S2 


20±5 


25 


Nsum 


SI 


742±64 


9 


R4k 


S2 


27±5 


18 


S6h 


S2 


19dz5 


25 










R41 


S2 


35±8 


22 


S6sum 


S2 


210±17 


8 











R4m S2 29±6 20 

masses in 10 b Mq. 



been calculated, the masses of the ionizing star clusters, 
M^ on , have been derived using the solar metallicity sin- 
gle burst models by Garcia- Vargas et al. (1995) which 
provide the number of ionizing photons per unit mass, 
[Q(Ho)/Mion]- A Salpeter initial mass function (Salpeter 
1955, IMF) has been assumed with lower and upper mass 
limits of 0.8 and 120 M©. In order to take into account the 
evolution of the Hll region, we have made use of the fact 
that a relation exists between the degree of evolution of the 



cluster, as represented by the equivalent width of the Hf3 
emission line, and the number of Lyman continuum pho- 
tons per unit solar mass (e.g. Diaz et al. 2000). We have 
measured the EW(H/3) from our spectra (see Table 6) fol- 
lowing the same procedure as in Hagele et al. (2006, 2008), 
that is defining a pseudo-continuum to take into account 
the absorption from the underlying stellar population. This 
procedure in fact may underestimate the value of the equiv- 
alent width, since it includes the contribution to the contin- 
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Figure 10. (cont) Enlargements of the F606W image around the 
CNSFRs of our study with the contours overlapped. The circles 
correspond to the adopted radius for each region. [See the elec- 
tronic edition of the Journal for a colour version of this figure 
where the adopted radii are in blue and the contours correspond- 
ing to the half light brightness are in red.] 



uum by the older stellar population (see discussions in Diaz 
et al. 2007 and Dors et al. 2008). The derived masses for the 
ionizing populations of the observed CNSFRs are given in 
column 6 of Table 6 and are between 1 and 7 per cent of the 
dynamical mass (see column 9 of the table). 

The amount of ionized gas (Mhii) associated to each 
star-forming region complex can also be derived from the 
Ha luminosities using the electron density (N e ) dependency 
relation given by Macchetto et al. (1990) for an electron 
temperature of 10 4 K. The electron density for each region 
(obtained from the [Sn] AA 6717 / 6731 A line ratio) has been 
taken from Pastoriza et al. (1993) for the regions in common. 
For the rest of the regions we assume a value of N e equal to 
100cm -3 , the value given by Pastoriza et al. (1993) for the 
other CNSFRs of this galaxy studied by them. For the nu- 
cleus of NGC 3310 we find an anomalous behaviour when we 
estimate Mhii using the values of the density (8200 cm -3 ) 
given by Pastoriza et al. (1993). We derive a rather low value 
of Mhii (5 x 1O 3 M0) for the Ha luminosity calculated by 
these authors. Then, when reviewing the density reported 
by Pastoriza and collaborators we found that it is consistent 
with the ratio of about 0.5 between the intensity values of 
the [Sll] AA 6717,6731 A emission lines given by them. How- 
ever, a careful inspection of the spectrum of the nucleus of 
NGC 3310 plotted in Fig. 5a of Pastoriza et al. (1993) reveals 
that the intensities of these two [Sn] emission lines are very 
similar to each other (ratio of about 0.9). A rough estimate 
using this revised ratio yields an electron density value of 
about 800 cm -3 . The corresponding Mhii value derived us- 
ing this lower density is an order of magnitude higher. 



6 DISCUSSION 

6.1 Star and gas kinematics 

Fig. 11 shows the relations between the velocity dispersions 
of gas and stars for NGC 3310. In the upper and lower panels 
we plot the relation derived from the Hf3 and [Oiii] emission 
lines, respectively. In both panels the straight line shows 
the one-to-one relation. In the upper panel the violet 3 cir- 
cles show the gas velocity dispersion measured from the Hf3 
emission line using a single Gaussian fit, maroon squares 
and blue downward triangles show the values measured from 
the broad and narrow components respectively using two- 
component Gaussian fits. The deviant points, marked with 
arrows, correspond to region R1+R2, which has a spectrum 
with low signal-to-noise ratio for which the cross-correlation 
method does not provide accurate results. In the lower panel, 
violet upward triangles, maroon diamonds and blue left tri- 
angles correspond to the values obtained by a single Gaus- 
sian fit, and to the broad and narrow components of the 
two Gaussian fits, respectively. Again, the deviant points 
(R1+R2) are marked with arrows. 

The H(3 velocity dispersions of the CNSFRs of 
NGC 3310 derived by a single-component Gaussian fit are 
the same, within the errors, as the stellar ones, except for 
R1+R2 and R7, for which they are lower by about 50 and 
25kms -1 , respectively. For R1+R2 this can be due to an 

3 In all figures, colours can be seen in the electronic version of 
the paper. 
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Table 6. Physical parameters. 



Region c(Ha) L(Ha) Q(H ) EW(H/3) M lon M H n M, on /M* 

(per cent) 

R1+R2 0.23 a 102.0 a 74~9 28i3 13^9 lOO 7 3738 L5 

R4 0.23 a 144.0 a 106.0 32.4 17.6 100 b 4.78 1.6 

R4+R5 0.20 a 218. a 160.0 41.7 21.4 100 b 7.24 1.1 

R6 0.17 a 57.3 a 42.1 16.7 12.4 100 c 1.90 3.5 

S6 0.35 b 62.5 b 45.9 12.5 17.4 100 b 2.07 6.6 

R7 0.17 a 45.5 a 33.5 19.4 8.7 100 c 1.51 1.0 

RIO 0.23 a 45. 5 a 33.5 9.7 15.7 100 c 1.51 2.4 

N 0.42 b 113.0 b 82.9 11.0 34.9 800 d 0.47 4.7 

J 0.23 a 573. a 421.0 82.5 31.4 200 b 9.52 

0.19 b 236.0 b 174.0 82.5 12.9 200 b 3.93 



Note. Luminosities in 10 38 ergs 1 , masses in 10 5 Mq, ionizing photons in 10 50 photons 1 
and densities in cm -3 
a From Diaz et al. (2000). 
b From Pastoriza et al. (1993). 

c Assumed from the values given by Pastoriza et al. (1993) for the other CNSFRs of this 
galaxy studied by them. 

^Derived from the spectrum of the nucleus of NGC 3310 plotted in Fig. 5a of Pastoriza 
et al. (1993) (see discussion in the text). 



overestimate of a*. Given the relatively low metal abun- 
dance of NGC 3310 (Pastoriza et al. 1993) the emission 
lines have generally a very good signal-to-noise ratio in the 
spectra of the CNSFRs, while for the red continuum it de- 
pends on the particular case (see Figs. 3 and 5 for the spa- 
tial profiles and the spectra, respectively). On the other 
hand, the stellar velocity dispersions of the CNSFRs are 
lower than those of the broad component of H(3 by about 
20kms -1 , again except for R1+R2 where a* is greater by 
about 35kms _1 , and for R7 where a* and H(3 broad are 
in good agreement. The narrow component of the CNSFRs 
shows velocity dispersions very similar for all the regions, 
with an average value of 23.0 =b 1.3 km s _1 , and is repre- 
sented as a blue solid line in the upper panel of Fig. 11, 
while the dashed lines of the same colour represent its error. 

The [Oiii] emission lines show the same behaviour as the 
H/3 lines for the two different components. In all cases the 
gas narrow component has velocity dispersion lower than 
the stellar one, and the broad component is slightly over 
it. The broad component of the [Oiii]5007A emission line 
(maroon diamonds in the figure) is wider than the stellar 
lines by about 20kms _1 , except for R4+R5 (S2 slit) and 
R7 for which they are approximately similar. The narrow 
component of [Oiii] (blue left triangles in the figure) shows a 
relatively constant value with an average of 20.9 ± 1.8 km s _1 
(blue lines in the figure). 

In general, the broad and narrow components of the H/3 
line have comparable fluxes. The ratios between the fluxes 
in the narrow and broad components of the H(3 emission 
line of the regions (including J, X and Y) vary from 0.95 to 
1.65, except for R10 and R7 for which we find ratios of 0.65 
and 2.5, respectively. The ratio of the narrow to broad [Oiii] 
fluxes is between 0.85 and 1.5; R10 and R7, for which this 
value is 0.65 and 1.98 respectively, are again the exception. 

The behaviour of the velocity dispersions of the CNS- 
FRs is found to be different in NGC 3310 from what was 
found for NGC 2903 and NGC 3351 (Papers I and II). For 
NGC 3310 we find a very good agreement in most cases be- 
tween velocity dispersions derived from gaseous and stellar 
lines assuming a single component for the gas. For NGC 2903 



and NGC 3351, in contrast, the velocity dispersions derived 
from emission lines measured using a single Gaussian fit 
are very different from the stellar ones (lower for H(3 and 
higher for [Oiii]) by about 20kms _1 ; see discussion in Pa- 
per II). Remarkably, the velocity dispersions derived from 
the narrow components of the two Gaussian fits show a rel- 
atively constant value for the whole CNSFR sample (about 
23kms- 1 ). If this narrow component is identified with ion- 
ized gas in a rotating disc, therefore supported by rotation, 
then the broad component could, in principle, correspond to 
the gas response to the gravitational potential of the stellar 
cluster, supported by dynamical pressure, explaining the co- 
incidence with the stellar velocity dispersion in most cases in 
the H(3 line (see Pizzella et al. 2004, and references therein). 
The velocity excess shown by the broad component of H(3 
and [Oiii] in most CNSFRs of NGC 3310 could be identified 
with peculiar velocities in the ionized gas (mainly in the 
high ionization gas) related to massive star winds or even 
supernova remnants. 

Our interpretation of the emission line structures, in the 
present study and in Papers I and II, parallels that of the 
studies of Westmoquette and collaborators (see for exam- 
ple Westmoquette et al. 2007a,b). They observed a narrow 
(~ 35-100 km s" 1 ) and a broad (~ 100-400 km s" 1 ) compo- 
nent to the Ha line across their four fields in the dwarf 
galaxy NGC 1569. They conclude that "the most likely ex- 
planation of the narrow component is that it represents 
the general disturbed optically emitting ionized interstellar 
medium (ISM), arising through a convolution of the stirring 
effects of the starburst and gravitational virial motions". 
They also conclude that "the broad component results from 
the highly turbulent velocity field associated with the inter- 
action of the hot phase of the ISM (material that is photo- 
evaporated or thermally evaporated through the action of 
the strong ambient radiation field, or mechanically ablated 
by the impact of fast-flowing cluster winds) with cooler gas 
knots, setting up turbulent mixing layers (e.g. Begelman & 
Fabian 1990; Slavin et al. 1993)". However, our broad com- 
ponent velocity dispersion values derived from H/3, all sig- 
nificantly lower than lOOkms -1 and similar to the stellar 
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velocity dispersions, resemble more their narrow component 
values. This is the same behaviour that we find in the broad 
component of the [Om] emission line for the CNSFRs of 
NGC 3310. The [Om] line profiles of the regions in Papers I 
and II on the other hand, behave as those described by West- 
moquette and collaborators. 

There are other studies that identified an underlying 
broad component to the recombination emission lines such 
as Diaz et al. (1987); Terlevich et al. (1996) in the M33 
giant Hn region NGC 604; Chu & Kennicutt Jr. (1994); 
Melnick, Tenorio-Tagle k Terlevich (1999) in the central 
region of 30 Doradus; Mendez & Esteban (1997) in four 
Wolf-Rayet galaxies, and Homeier & Gallagher (1999) in the 
starburst galaxy NGC 7673. More recently, Westmoquette et 
al. (2007c, 2009) found a broad feature in the Ha emission 
lines in the starburst core of M82; Ostlin et al. (2007) and 
James et al. (2009) in the blue compact galaxies ESO 338- 
IG04 and Mrk996, respectively; Sidoli et al. (2006); Firpo 
et al. (2009) in giant extragalactic Hn regions, and Hagele 
(2008) in CNSFRs of early type spiral galaxies. The first two 
and the last three studies also found this broad component 
in the forbidden emission lines. 

In the nuclear region of the galaxy the stars and the 
broad component of the ionized gas, both for Hf3 and [Om] 
have the same width, with a a of about 73 and 83kms _1 , 
respectively. The profiles of the nuclear lines hence, behave 
similarly to what is described by Westmoquette and collabo- 
rators. The gas narrow component however shows a substan- 
tially lower velocity dispersion though still higher by about 
15kms~\ than that of the CNSFRs. 

For regions J, X and Y we can not derive a stellar ve- 
locity dispersion due to the low signal-to-noise ratio of their 
continuum and the noise added by subtracting the large 
amount of emission lines present in their red spectra (see 
right panels of Fig. 5). However, if we analyze the velocity 
dispersions derived from their strong emission lines in the 
blue spectral range we find values very similar to those of 
the other CNSFRs (see Table 2). 

Fig. 12 presents the radial velocities derived from the 
H/3 and [Om] emission lines and the CaT absorptions along 
the slit for each angular slit position of NGC 3310, SI in the 
upper panel and S2 in the lower one. The rotation curves 
seem to have the turnover points located in or near the 
star-forming ring. Due to the relatively low metallicity of 
NGC 3310 (Pastoriza et al. 1993), the gas emission lines 
are very strong in the central zone of the galaxy and we 
can derive the radial velocity of the gas much further (up 
to 18") than for the stellar CaT absorptions. The gas 11(3 
and [Om] 5007 A radial velocities and the stellar ones are in 
very good agreement, except in the zone located around R4, 
where we find a difference between gas and stars of about 30 
to 40kms _1 . These differences could be due to a low signal- 
to-noise ratio of the stellar continuum emission in this zone 
of the galaxy. 

In both slit positions we can appreciate a very steep 
slope in the radial velocity curve as measured from H/3 and 
[Om], in very good agreement with each other. These strong 
gradients occur more or less at the same distance from the 
nucleus (~8") towards the north-east of the galaxy, almost 
at the position where the two slits intersect (see Figs. 1 and 
2). The curves show a step in the radial velocity of about 
70 and 90kms _1 for SI and S2, respectively. In Fig. 2 we 



indicate with two circles of 8 and 18" radii the approximate 
position where this step is located and the distance from the 
nucleus up to where we can derive the radial velocities for 
SI. The stellar radial velocity of R1+R2 derived using the 
5 px aperture (criss-crossed white-blue diamond in the lower 
panel of Fig. 12) is in very good agreement with the values 
derived from the gas emission lines. 

A very similar result was found by Pastorini et al. (2007) 
using high spatial resolution spectra (~ 0.07" px -1 ) and 
moderate spectral resolution (R= A A/ A ~ 6000, yielding a 
resolution of AA ~ 1.108 A px -1 at Ha) from the Space Tele- 
scope Imaging Spectrograph (STIS) on board the HST, and 
using a 2 x 2 on-chip binning for the non nuclear spectra. 
The position angle of their slits is 170° and three parallel 
positions were used, one of them passing across the nucleus 
and the other two with 0.2 " offset to both sides 4 . As pointed 
out by Pastorini and collaborators, NGC 3310 shows a typ- 
ical rotation curve expected for a rotating disc (a typical S 
feature, Marconi et al. 2006) but the curve is disturbed. Be- 
sides, we can appreciate that near the Jumbo region there 
is a deviation from circular rotation motion. The step in the 
radial velocity found by us is an effect of the spatial resolu- 
tion. In Fig. 3 of Pastorini et al. (2007), the better spatial 
resolution of STIS-HST resolves the steep gradient and com- 
pletes the information missing in our data. This disturbed 
behaviour of the rotation curve of the gas in the galactic 
disc, characterized by non-circular motions, was also found 
by Mulder et al. (1995) using Hi radio data; Mulder & van 
Driel (1996) from Hi and Hn radio data and Kregel & San- 
cisi (2001) using optical data. They found a strong streaming 
along the spiral arms, supporting the hypothesis of a recent 
collision with a dwarf galaxy that triggered the circumnu- 
clear star formation during the last 10 8 yr. Besides, van der 
Kruit (1976) shows that the rotation centre of the gas is off- 
set with respect to the stellar continuum (by ~ 1.5"). The 
central velocity of NGC 3310 derived by us is in very good 
agreement with that previously estimated by van der Kruit 
(1976); Mulder et al. (1995); Mulder & van Driel (1996); 
Haynes et al. (1998); Kregel & Sancisi (2001). The velocity 
distribution is also in very good agreement with that ex- 
pected for this type of galaxies (Binney & Tremaine 1987). 

Fig. 13 shows the run of velocity dispersion along the 
slit versus pixel number for slit positions SI and S2 respec- 
tively. These velocity dispersions have been derived from 
the gaseous emission lines, H/3 and [Om], and the stellar 
absorption ones using 2 and 3px apertures for SI and S2, 
respectively. We have marked the location of the studied 
CNSFR and the galaxy nucleus. We have also plotted the 
stellar velocity dispersion derived for each region and the 
nucleus using 5 px apertures. The values derived with wider 
apertures are approximately the average of the velocity dis- 
persions estimated using the narrower ones. Therefore, the 
apparent increase of the velocity dispersions of the regions 



4 It must be noted that the Ha narrow band-filter image in Fig. 
2 of that work (the same ACS image presented by us) is rotated 
by 180° with respect to their quoted position, since the Jumbo 
region appears placed at the north-east of the nucleus (according 
to the orientation of the image given by these authors) while this 
region is located at the south-west (see e.g. Telesco Sz Gatley 
1984; Terlevich et al. 1990; Pastoriza et al. 1993). 
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lying on relatively steep parts of the rotation curve by spa- 
tial "beam smearing" of the rotational velocity gradient due 
to the finite angular resolution of the spectra is not so criti- 
cal for the studied CNSFRs. This is due to the fact that the 
turnover points of the rotation curves and the star-forming 
ring seem to be located at the same positions, as was pointed 
out above. The velocity dispersions, both for stars and gas, 
show a behaviour characteristic of a regular circular motion 
in a rotating disc, a smooth plateau and a central peak. A 
similar result was found by Pastorini et al. (2007). 

6.2 Star cluster masses 

The eight CNSFRs observed in NGC 3310 show a complex 
structure at the HST resolution (Fig. 10), with a good num- 
ber of subclusters with linear diameters between 2.2 and 6.2 
pc. For these individual clusters, except for R1+R2, the de- 
rived upper limits to the masses are in the range between 
1.8 and 7.1 x 10 6 M0 (see table 5), with fractional errors be- 
tween rsj 16 and ~ 32 per cent. The upper limits to the dy- 
namical masses estimated for the whole CNSFRs ("sum"), 
except for R1+R2, are between 2.1 x 10 7 and 1.4 x 10 8 M , 
with fractional errors between ~3 and ^8 per cent. The 
stellar velocity dispersion of region R1+R2 has a large asso- 
ciated error, so we do not quote an explicit error for its de- 
rived mass which we consider highly uncertain. The masses 
for the individual clusters of this CNSFR, derived using this 
value of the velocity dispersion are larger, between 1.2 and 
2.4 x 1O 7 M0. However, if we assume that the relation be- 
tween the gas and the stellar velocity dispersions follows a 
behaviour similar to that shown by the rest of the CNSFRs 
of NGC 3310 (see Fig. 11) these masses could be reduced by 
a factor of 2.3, yielding individual masses between 5.2 x 10 6 
and 1.0 x 10 7 M o . In the same way, the total dynamical mass 
of R1+R2, 9.1 x 10 8 M , would be reduced to 4.0 x 1O 8 M . 
The upper limits to the dynamical mass derived for the nu- 
clear region inside the inner 14.2 pc is 5.3 x 1O 7 M , with a 
fractional error of about 11 per cent. This value is some- 
what higher than that derived by Pastorini et al. (2007) 
under the assumption of the presence of a super massive 
black hole (SMBH) in the centre of the galaxy. When they 
allow for the disc inclination, this mass varies in the range 

5.0 x 10 6 -4.2x 1O 7 M . 

The masses of the ionizing stellar clusters of the CNS- 
FRs, derived from their Ha luminosities (under the assump- 
tion that the regions are ionization bound and without tak- 
ing into account any photon absorption by dust), are be- 
tween 8.7 x 10 5 and 2.1 x 10 6 M for the star-forming re- 
gions, and 3.5 x 10 6 M for the nucleus (see Table 6). In the 
case of the Jumbo region this mass vary between 1.3 and 

3.1 x 10 6 M whether we take the Ha luminosities from Pas- 
toriza et al. (1993) or Diaz et al. (2000) respectively, proba- 
bly because the estimated reddening constants in these two 
works are different (see Table 6). These different values of 
the luminosities can be due to the very complex structure 
of the region (see Fig. 1 of both works) and to different se- 
lection of the zone used to measure the Ha fluxes. Besides, 
Pastoriza and colleagues used spectroscopic data while Diaz 
and collaborators used photometric images. In column 9 of 
Table 6 we show a comparison (in percentage) between the 
ionizing stellar masses of the circumnuclear regions and their 
dynamical masses. These values are approximately between 



1-7 per cent for the CNSFRs, and 5 per cent for the nucleus 
of NGC 3310. 

Finally, the masses of the ionized gas, also derived from 
their Ha luminosities, range between 1.5 and 7.2 x 10 5 M 
for the CNSFRs, and 4.7 x 10 4 M for the nucleus (see Table 
6). They make up a small fraction of the total mass of the 
regions. As in the case of the masses of the ionizing stellar 
clusters of the Jumbo region and for the same reasons, the 
mass of the ionized gas derived using the Ha luminosities 
from Pastoriza et al. (1993) or from Diaz et al. (2000) are 
different, 3.93 and 9.52 x 1O 5 M , respectively. It should be 
taken into account that we have derived both the masses of 
the ionizing stellar clusters and of the ionized gas from the 
Ha luminosity of the CNSFRs assuming that they consist 
of one single component. However, if we consider only the 
broad component whose kinematics follows that of the stars 
in the regions, all derived quantities would be smaller by a 
factor of 2. 

Comparing the colours of the small-scale clusters from 
HST data with models, Elmegreen et al. (2002) suggest that 
most of these clusters contain masses from 10 4 M to sev- 
eral times 1O 4 M . For the large-scale "hot-spots" they es- 
timate masses ranging from 10 4 to several times 1O 5 M . 
For their 108+109 region (almost equivalent to the Jumbo 
region) they derived a mass of 6.3 x 10 5 M , similar to the 
value derived by Pastoriza et al. (1993), 7x 1O 5 M . They 
estimated from the unusually large H and Hell emission line 
luminosities, that this region must contain 220 WN and 570 
OB stars. Regarding other regions, Elmegreen and colabo- 
rators derived a mass of 10 6 M for their region 78 using 
the NIR data and 7.7 x 10 6 M from their optical ones. Be- 
sides, they found 17 candidate SSCs, with masses in the 
range 2 x 10 4 to 4.5 x 10 5 M©. Some of them coincide with 
the CNSFRs studied by us, such as 48, 49 and 50, approxi- 
mately coincident with R5 and R6. These SSCs are mostly in 
the innermost southern spiral arm, with some in the north- 
ern one or outside the southern arm of the ring (Elmegreen 
et al. 2002). Their derived (J-H) colours suggest two differ- 
ent populations of SSCs, a very young of few million years 
and an older one in a range between 10 and 50 million years, 
with the younger clusters located in the northern part of the 
ring. 

It should be recalled that we have estimated the dynam- 
ical masses through the virial theorem under the assumption 
that the systems are spherically symmetric, gravitationally 
bound and have isotropic velocity distribution. We have used 
the stellar velocity dispersions derived from the CaT absorp- 
tions features and the cluster sizes measured from the high 
spatial resolution WFPC2-PC1 HST image. As mentioned 
at the end of section 4, the use of these wider size scale 
velocity dispersion measurements to estimate the mass of 
each knot, can lead to an overestimate of the mass of the 
individual clusters, and hence of each CNSFR (see Paper I). 

However, as can be seen in the HST-NICMOS images 
(right-hand panel of Fig. 1), the CNSFRs are clearly visi- 
ble in the IR and dominate the light inside the apertures 
observed (see Elmegreen et al. 2002, for a detailed analy- 
sis of the IR images) . All the regions analyzed show up very 
prominently in the near-IR and therefore we can assume that 
the light at the CaT wavelength region is dominated by the 
stars in the clusters. The IR CaT is very strong, in fact the 
strongest stellar feature, in young clusters older than 4 Myr 
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Terlevich et al. (1990). Besides, we detect a minimum in the 
velocity dispersion at the position of most of the clusters, 
indicating that they are kinematically distinct. We can not 
be sure though that we are actually measuring their velocity 
dispersion and thus prefer to say that our measurements of 
a* and hence the dynamical masses constitute upper limits. 
Although we are well aware of the difficulties, still we are 
confident that these upper limits are valid and important 
for comparison with the gas kinematic measurements. 

Another important fact that can affect the estimated 
dynamical masses is the presence of binaries among the red 
supergiant and red giant populations from which we have de- 
rived the stellar velocity dispersions. In a recent work, Bosch 
et al. (2009) using the Gemini Multi-Object Spectrograph 
(GEMINI-GMOS) data, have investigated the presence of 
binary stars within the ionizing cluster of 30 Doradus. From 
a seven epoch observing campaign they have detected a rate 
of binary system candidates within their OB star sample 
of ^50%. Interestingly enough, this detection rate is con- 
sistent with a spectroscopic population of 100 % binaries, 
when the observational parameters described in the simula- 
tions by Bosch & Meza (2001) are set for their observations. 
From their final sample of 'single' stars they estimated a ra- 
dial velocity dispersion of 8.3 km s -1 . When they derived a* 
from a single epoch, they found values as high as 30kms _1 , 
consistent with the values derived from single epoch NTT 
observations by Bosch et al. (2001). 



Although the environment of our CNSFRs is very dif- 
ferent from that of 30 Dor and the stellar components of 
the binary systems studied by Bosch et al. (2009) are very 
different from the stars present in our regions from where 
the CaT arise (red supergiants) , this is an illustrative ob- 
servational example of the problem. The orbital motions of 
the stars in binary (multiple) systems produce an overesti- 
mate of the velocity dispersions and hence of the dynamical 
masses. The single-star assumption introduces a systematic 
error that depends on the properties of the star cluster and 
the binary population, with an important effect on the clus- 
ter mass if the typical orbital velocity of a binary component 
is of the order of, or larger than, the velocity dispersion of the 
single/binary stars in the potential of the cluster (Kouwen- 
hoven & de Grijs 2008). As was pointed out by these authors, 
the relative weights between the single and binary stars in 
the velocity dispersion measurements depend on the binary 
fraction, which, together with the semi-major axis or period 
distribution, are the most important parameters in order 
to determine if the binary population affects the estimated 
dynamical masses. Their simulations indicate that the dy- 
namical mass is overestimated by 70%, 50% and 5 % for a 
measured stellar velocity dispersion in the line of sight of 
lkms -1 , 2kms _1 and lOkrns -1 respectively. They there- 
fore conclude that most of the known dynamical masses of 
massive star clusters are only mildly affected by the presence 
of binaries. Hence, although the binary fraction of the red 
supergiants and red giants in our circumnuclear regions is 
unknown, since the smallest estimated velocity dispersion is 
31kms _1 , we can assume that the contribution of binaries 
to the stellar velocity dispersions is not important. 



7 SUMMARY AND CONCLUSIONS 

We have measured gas and stellar velocity dispersions in 
eight CNSFRs and the nucleus of the barred spiral galaxy 
NGC 3310. The stellar velocity dispersions have been mea- 
sured from the CaT lines at AA8494, 8542, 8662 A, while 
the gas velocity dispersions have been measured by Gaus- 
sian fits to the A 4861 A and the [Om] A5007A emission 
lines on high dispersion spectra. 

Stellar velocity dispersions are between 31 and 
73kms _1 . The stellar and gas velocity dispersion are in 
relatively good agreement, with the former being slightly 
larger. The velocity dispersions from [Oiii]5007A behave 
similarly to those from H/3. However, the best Gaussian fits 
involve two different components for the gas: a "broad com- 
ponent" with a velocity dispersion larger than that mea- 
sured for the stars by about 20kms _1 , and a "narrow com- 
ponent" with a velocity dispersion lower than the stellar one 
by about 30kms _1 . This last velocity component shows a 
relatively constant value for the two gas emission lines, close 
to 22kms _1 and also close to that measured for CNSFRs in 
the previously studied NGC 3351 and NGC 2903. The ve- 
locity dispersion for the broad component of H/3 is 
similar to that of [OIII], contrary to what we found 
in Papers I and II for NGC 2903 and NGC 3351. We 
find a velocity shift between the narrow and broad compo- 
nents of the multi-Gaussian fits that vary between -25 and 
20kms _1 in radial velocity. 

When plotted in a [Olll]/H/3 versus [Nil] /Ha diagnostic 
diagram, the broad and narrow systems and those values 
derived using a single Gaussian fit show very similar values, 
lying in the region of low-metallicity Hn-like objects. This 
is in contrast to what we found in Papers I and II where the 
two systems are clearly segregated for the high-metallicity 
regions of NGC 3351 and NGC 2903, with the narrow com- 
ponent showing lower excitation and being among the lowest 
excitation line ratios detected within the SDSS data set of 
starburst systems. 

The rotation curve of NGC 3310 shows a typical S 
feature, with the presence of some perturbations, partic- 
ularly near the location of the Jumbo region. The values 
derived from the gas emission lines and the stellar absorp- 
tion features are in very good agreement. The position going 
through the nucleus shows maximum and minimum values 
more or less coincident with the location of the CNSFRs. 

The upper limits to the dynamical masses estimated 
from the stellar velocity dispersion using the virial theo- 
rem for the CNSFRs are in the range between 2.1 x 10 7 and 
1.4 x 10 8 Mq. For the nuclear region inside the inner 14.2 pc 
this upper limit is 5.3 x 10 7 M©. The upper limits to the de- 
rived masses for the individual clusters are between 1.8 and 
7.1 x 10 6 M . 

Masses of the ionizing stellar clusters of the CNSFRs 
have been derived from their Ha luminosities under the as- 
sumption that the regions are ionization bound and without 
taking into account any photon absorption by dust. These 
masses derived for the star forming complexes of NGC 3310 
are between 8.7 x 10 5 and 2.1 x 10 6 M , and is 3.5 x 10 6 for 
the nucleus (see table 6). Therefore, the ratio of the ionizing 
stellar population to the total dynamical mass is between 
0.01 and 0.07. 

Derived masses for the ionized gas, also from their Ha 
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luminosities, vary between 1.5 and 7.2 x 10 5 M© for the 
CNSFRs. For the nucleus, the derived mass of ionized gas is 
4.7 x 1O 4 M . 

It is interesting to note that, according to our findings, 
the SSC in CNSFRs seem to contain composite stellar pop- 
ulations. Although the youngest one dominates the UV light 
and is responsible for the gas ionization, it constitutes only 
about 10 per cent of the total mass. This can explain the low 
EWs of emission lines measured in these regions. This may 
well apply to studies of other SSC and therefore conclusions 
drawn from fits of single stellar population models should 
be taken with caution (e.g. McCrady et al. 2003; Larsen 
et al. 2004). Furthermore, the composite nature of the CNS- 
FRs means that star formation in the rings is a process that 
has taken place over time periods much longer than those 
implied by the properties of the ionized gas. 
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Figure 11. Upper panel: relation between velocity dispersions of 
the gas (derived from H/3) and stars (CaT) for the CNSFRs and 
the nucleus of NGC 3310. Symbols are as follows: Upper panel: 
single Gaussian fit, violet circles; two Gaussian fit, broad compo- 
nent, maroon squares; narrow component, blue downward trian- 
gles; lower panel: as the upper panel for the [Olll] line. Violet up- 
ward triangles correspond to the estimates using a single Gaussian 
fit, maroon diamonds represent the broad components of the two 
Gaussian fit and blue left triangles, the narrow components. The 
blue solid line represents the average velocity dispersion of the 
narrow component of the gas (Hf3 upper and [Olll] lower panel), 
for the CNSFRs, and the blue dashed lines represent their esti- 
mated errors. [See the electronic edition of the Journal for a colour 
version of this figure.] 
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Figure 12. Radial velocities along the slit versus pixel number 
for each slit position of NGC3310 (upper panel: SI; lower panel: 
S2) as derived from the gas emission lines (red circles: H/3; upward 
triangles: [Om]) and the stellar absorption ones (blue squares). 
The criss-crossed white-blue diamond in the lower panel corre- 
spond to the stellar velocity of R1+R2 derived using the 5px 
aperture. The individual CNSFRs and the nucleus, "N", or the 
closest position to it, "d", are marked in the plots. The dashed- 
dotted line is the velocity of the nucleus of NGC 3310 derived by 
van der Kruit (1976). The distance in arcsec from the nucleus is 
displayed in the upper x-axis of each panel. [See the electronic 
edition of the Journal for a colour version of this figure.] 
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Figure 13. Velocity dispersions along the slit versus pixel num- 
ber for each slit position of NGC 3310 (upper panel: SI; lower 
panel: S2) as derived from the gas emission lines (small red cir- 
cles: H/3; small yellow triangles: [Om]) and the stellar absorption 
ones (small blue squares). The stellar velocity dispersions derived 
for each region and the nucleus using the 5px aperture are also 
plotted with violet squares. The individual CNSFRs and the nu- 
cleus, "N", or the closest position to it, "d", are marked in the 
plots. The distance in arcsec from the nucleus is displayed in the 
upper x-axis of each panel. [See the electronic edition of the Jour- 
nal for a colour version of this figure.] 



